Charged Snowball in Non-Polar Liquid 
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The problem of correct definition of the charge carrier effective mass in superfiuid helium is 
revised. It is demonstrated that the effective mass M of a such quasiparticle can be introduced 
without use of the Atkins's idea concerning the solidification of liquid He in the close vicinity of 
ion 0. The two-liquid scenario of the "snowball" mass formation is proposed. The normal fluid 
contribution to the total snowball effective mass, the physical reasons of its singularity and the way 
of corresponding regularization procedure are discussed. Agreement of the theory with the available 
experimental data is found in the wide range of temperatures. 

PACS numbers: 67.40.-w, 72.20.Jv 



The ion-dipole interaction between the inserted 
charged particle and the induced electric dipoles of sur- 
rounded atoms is one of the interesting phenomena which 
take place in a non-polar liquid. It is this interaction 
which is generally responsible for the different salvations 
phenomena there |2j, while in non-polar cryogenics liq- 
uids (like He, Ne, Ar, etc.) it leads to so-called "snow- 
ball effect". The latter consists in the formation of a 
non-uniformity Sp (r) in a density of a liquid around the 
inserted charged particle. 

The ion-dipole interaction Ui-d (r) between the in- 
serted charged particle and solvent atoms in the simplest 
form can be written as 



Ui-d (r) 



a e 



(1) 



where a is the polarization of a solvent atom and r is the 
distance to the charged particle. The presence of attrac- 
tion potential in equilibrium has to be compensated 
by the growth of the solvent density Sp (r) in direction of 
the charged particle. The latter can be estimated, using 
the requirement of the chemical potential constancy: 



Sp (r) = P s (r) v. 

■dp_ 
dP 
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where P s (r) is the local pressure around the ion, v s is 
the volume of individual solvent atom and s is the sound 
velocity. The distance at which P s (r) reaches the value 
of pressure of a solvent solidification P/ 



R„ = 



2 P;v„ 



(3) 



corresponds to the radius of rigid sphere which is called 
"snowball". In the case of positive ions being in liquid 
helium, where the inter-atomic distance a«34 and the 



pressure P s s as 25 atm, the snowball radius is estimated 
as Pf e w 7 A. 

Described above so-called Atkins's snowball model 
is quite transparent and it was found useful for the var- 
ious qualitative predictions. In particular, it provides 
by the natural definition and estimation for the effective 
mass M of a such quasi-particle as the sum of the extra 
mass caused by the presence of the density perturbation 
Sp(r) 



M (soi) = 4?r / r 2s p (r)dr 



(4) 



and associated mass, related to the appearance of the ve- 
locity distribution around the sphere moving in an ideal 
liquid 



A^ ass) (p,R s ) = ^pR 3 s 



(5) 



Both of these contributions turn out to be of the same 
order. Calculated in this way snowball effective mass 
turns out to be M > 50m4, what roughly corresponds to 
the experimental data 

Careful analysis shows that both expressions and 
(0 require more precise definitions and further develop- 
ment of the Atkins's model. Indeed, one can easily see 
that the value of M^ so1 ^ turns out to be critically sensitive 
to the lower limit of the integral. Atkins used the value 
of helium inter-atomic distance a sa 3 ^4 .1} as a rough 
cut-off parameter only. The microscopic analysis shows 
(see Eq. (I2(J|) below) that the real lower limit turns out 
to be less that this value. 

Principal revision requires the definition of the asso- 
ciated mass M( ass ) for the two-fluid model. The matter 
of fact that its strong temperature dependence, which 
has been systematically observed experimentally 0, 
cannot be explained within the "ideal" flow picture. It is 
why below we propose the hydrodynamics scenario of the 
associated mass M^ ass ^ formation that takes into account 
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the viscous part of this problem. The straightforward ac- 
counting for the non-zero viscosity results in a dramatic 
growth (compared to the ideal case) of kinetic energy of 
the moving sphere and, consequently, its effective asso- 
ciated mass. This fact is caused by setting in motion of 
a spacious domain of viscous liquid around the moving 
particle. Nevertheless, the arising divergency can be cut 
off by accounting for the non-linear effects in stationary 
flow. 

Another way to introduce the effective associated mass 
for the problem under consideration is based on the anal- 
ysis of the real experiments where the dynamical pa- 
rameters of charged carriers are always measured in os- 
cillating regime. In this approach the problem can still be 
assumed linear but the effective mass becomes frequency 
dependent. The proposed consideration permits us to 
explain consistently the specific features of the observed 
on experiments snowball effective mass temperature de- 
pendence in wide interval of temperatures. 

Let us start from the definition of the effective asso- 
ciated mass within the "snow-cloud" approximation and 
using the two-fluid model. In that case 



M (ass) = M ( aSS ) +M) (ass)_ 

The first term can be defined as follows: 



s)y2 /2= P^Zl /°°v2(r)d 3 r, 
2 Jr. 
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where p s (T) is the super-fluid density, V is the velocity 
of a "snow-cloud" center of mass forward motion, v s (r) 
is the super-fluid component local velocity distributions 
appearing in liquid due to charge carrier motion. Let 
us underline that some auxiliary parameter i?* instead 
of the snowball radius R^ c (see Eq.@) here is used as 
the lower limit cut-off. Its origin and evaluation will be 
discussed below (see Eq.(j20J)). 

The super-fluid flow has a potential character d, Q 



Vcj> (r) , (r) 



A n 



A = VRl/2. (8) 



Correspondingly the super-fluid part of associated mass 
Ms ass ^ with the velocity distribution determined by 
the Eq. JSJ is reduced to the expression JSJ with the 
simple substitutions R s — > R* and p — > p s . 

As to the normal fluid part of the associated mass 
Mn ass \ its definition turns out to be more cumbersome. 
It would be natural to define it in the same way as the 
superfluid one, just substituting in Eq. (0) the subscript 
s by n. But this programme runs against the consid- 
erable obstacles in spite of the fact that the problem of 
rigid sphere motion in a viscous liquid was considered 
a long ago (see, for instance, [1,0). Indeed, for small 



Reynolds numbers Re = pVR^/r] <C 1 and distances 
-R* < r < i?»/Re the gradient term (v • V) v in cor- 
responding Navier-Stokes equation can be ignored. The 
solution of such linearized equation in zero approxima- 
tion turns out to be independent on the liquid viscosity 
j]. In the center of mass frame it has the form: 



(9) 
(10) 



Here 6 is the polar angle counted from the x-axes, which 
coincides with the ion velocity direction. 

One can easily see, that corresponding contribution to 
the kinetic energy diverges at the upper limit of integra- 
tion. To regularize this divergence it is necessary to use at 
large distances r>l r) = R*/ Re (l n = rj/ pV is the charac- 
teristic viscous length) more precise, so-called Ossen, so- 
lution of the Navier-Stokes equation which was obtained 
taking into account the gradient term in it [9j . This solu- 
tion shows that almost for all angles, besides the domain 
restricted by the narrow paraboloid (x) = ir—^Jrj/pV\x\ 
behind the snowball (so-called "laminar trace") at large 
distances the velocity field decays exponentially. In latter 
the velocity decays by power law and it gives logarith- 
mically large contribution with respect to other domain 
of disturbed viscous liquid. Such specification permits 
to cut off formal divergency of the kinetic energy and 
to find with the logarithmic accuracy the value of the 
normal component of associated mass for the stationary 
moving in viscous liquid charge carrier: 
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Here the associated mass Mq 8 ' (p„,i?*) is defined by 
Eq. JSJ) with radius i?* and density p n . Since in our as- 
sumptions l v 3> i?* the value of M^f^l turns out much 
larger than the value of the associated mass in the ideal 
liquid. Moreover, when velocity V — > such definition, 
in spite of the performed above regularization procedure, 
fails since Ea. l|llJ) diverges. 

At this point it is necessary to recall that, as it was 
mentioned in introduction, in reality it is the dynamic 
associated mass which has to be considered. This cir- 
cumstance allows one to incorporate qualitatively new 
features in the effective mass behavior. Indeed, let us 
consider the situation, when a periodic electric field is 
applied to a charge carrier placed in liquid He 4 . The dy- 
namic Stokes force, appearing when the sphere oscillates 
in the viscous liquid with finite frequency, has the form 
7, 8]: 



(11) 
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F (oj) = 67r?7i?„ 1 



R* 



+ 3^/^(1 



95 (w) 



(12) 



(to) , 



1/2 

where <5(w) = (2ri/p n uj) is so-called dynamic penetra- 
tion depth. It is natural to identify the coefficient in front 
of the Fourier transform of acceleration iojV (oj) with the 
effective dynamic associated mass, what gives: 

<t) H = ^ SS) (Pn,R.) (l + l 5 -^) • (13) 

One can see that for high frequencies (6 (oj) <C i?*) the 
dynamic associated mass coincides with that one of a 
sphere moving stationary in an ideal liquid, while when 



OM 



(ass) 



n(dyn 



(oj) diverges as lo 1 / 2 . As we already have 
seen above this formal divergence is related to fall down of 
the linear approximation in the Navier-Stokes equation 
assumed in derivation of Ea.Q12[l. It is clear that the 
definition l|13|) is valid for high enough frequencies, until 
5 (lo) < l v (i.e. oj >lu — ?//p„i?^).When oj becomes lower 
than oj the penetration depth 5 (lo) in Ea. i|13|) has to be 
substituted by l v and up to the accuracy of In l v / R* the 
dynamic definition Ea. (|13|) matches with the static one 
(see Eq.JTTJl). 

Thus, the following qualitative picture arises. In the 
stationary regime, the moving sphere disturbs an enor- 
mous domain in viscous liquid. This perturbation results 
in the large effective mass (|11|) and can substantially af- 
fect the critical conditions for the generation of vertex 
rings. When going to the oscillatory regime the viscous 
"coat" and especially the "laminar trace" do not keep 
pace with the charge in its oscillations, the perturbed 
domain shrinks with the growing frequency thus reduc- 
ing the value of effective mass. 

The question arises, how sensitive are the introduced 
above definitions of the associated mass Eas. (|13|) . Q) to 
the real shape of the liquid density perturbation in the 
vicinity of the charge carrier? The examples are already 
known when the snowball Atkins's model and more re- 
alistic snow-cloud model lead to qualitatively different 
predictions for the value of positive ion mobility in liq- 
uid He 4 (see 0-^3 )• I n order to clear up this problem 
let us consider the hydrodynamic picture of a snow-cloud 
motion. The latter we assume as the density compression 
which decays with distance by power law: 



p (r) = p + Sp (r) = p + 



C_ 

^4 



(14) 



( constant C is expressed in terms of Ecis.©). The con- 
tinuity equation is read as 



p (r) diw + v • Vp (r) = 0. 



(15) 



For an ideal (superfluid) liquid the velocity field is po- 
tential and it can be presented in the form Eq.JSJ). Sub- 
stituting Eqs.JTlI) and © to Eq.JTSJ one finds 



Acj)(r) 



AC (r) 
r 5 p + Sp 



0. 



(16) 



The Eq. (|16|) has to be solved with the additional require- 
ment 



' (r — > oo) — > rV cos ( 



(17) 



In Born approximation the Ea. (|16|l is reduced to the 
Poisson equation. Supposing that </> (r) — rV cos 6 + 
4>i (r) , it can be rewritten in the form 



AC V cos6 , 
r D p + op 

which solution can be presented as 



01 (x) 



d-: 



7T 



(Vt) 



? ,e p + Sp (r) 



(19) 



The Ea. (|19(l shows, that far enough from a snow-cloud 
center the potential of velocity distribution looks like a 
dipole one, exactly in the same way as for the flow around 
the rigid sphere of the radius i?idcai- Using Eqs.© and 
(fTU|l one finds 



I? 3 



deal 



2ttC 



dr 



r 2 (p + Sp) 



(20) 



It is easy to see that in spite of presence of r 2 in the 
denominator this integral converges at the lower limit. 

The above analysis demonstrates that the problem of 
the associated mass definition in the snowball and the 
snow-cloud models for an ideal liquid are qualitatively 
identical. It is just enough to renormalize the effective 
radius of the hydrodynamic dipole. Elsewhere will be 
shown that the same conclusion can be done basing on 
the Navier-Stokes formalism for the motion of the snow- 
cloud in the viscous liquid (naturally, with other defini- 
tion of the effective radius). 

Let us pass to discussion of the experimental situation. 
First of all it is necessary to stress that the effective mass 
temperature dependence for positive charge carriers in 
liquid He 4 has been observed in the interval from 50 mK 
up to 2 if (see 0-0 )• The growth of the mass with 
temperature occurs with the average rate 28 m He 4 / K. 
The viscous scenario does not permit us to analyze the 
region of very low temperatures, so let us restrict our 
consideration by the interval 0.5 K < T < 2 K, where 
the density of the normal component varies in the limits 
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0.005 < p n /p < 0.7. The data concerning temperature 
dependence M< ass > (T) were taken from the ac transport 
measurements. Namely, the authors of Ref. Q| analyzed 
temperature dependence of the position of maximum ob- 
served in frequency dependent conductivity. This posi- 
tion in our model evidently corresponds to situation when 
both the real and imaginary parts of the Stokes force l|12|) 
are of the same order. Equating them one can evaluate 
the corresponding frequency as 



w x (T) 



x 2 rj (T) 



9/2. 



2p n (T) Rl 

Corresponding value of the associated mass is 



AI 



(ass) 
(dyn 



K) = 3-KR^pn- I + t:^ 



(21) 



(22) 



(ass) 



Fitting of the Eq. © with the Eqs. Q [HJj for M, 
and ilL'21) for M n ass ^ demonstrates a good agreement with 
the experimentally observed temperature dependence in 
the wide range of temperatures (see Fig. 0. This fact 
strongly supports the idea to use two-liquid scenario and 
to account for the viscous velocity distribution in the 
frameworks of the Navier-Stokes flow picture. It is nec- 



120 



90 - 



60 - 



30 f 2 *¥ a a - - a - 

* A - . . 

,U • • ' . 



1,3 



1,5 



1.9 



T(K) 



FIG. 1: The temperature dependencies of Mi ass ' (T) accord- 
ing to Eqs. 0EJ (triangles), M^ Bs) (T) (see Eq. (22), black 
dots) and their sum Afj ass) (T) + Mi™ s) (T) (crosses) are pre- 
sented separately. The dependence p n (T) is taken from |14|. 
The experimental data (open circles with error bars) are taken 
from Ref. Q 

essary to mention that understanding of possible viscous 
origin of the temperature dependence (T) was al- 

ready contained in the e arly paper Q. Without deriva- 
tion, just referring to 0.1131 > the authors presented their 
definition of the associated mass, which resembles much 
our formula (| 1 31) . but contained the additional factor 
p n / P in its second term. This difference turns out es- 
sential for the fitting of experimental data, leading to 
the visible discrepancy at low temperatures. That is per- 
haps the reason why in Ref. |4J the real M{T) fitting is 



only realized in a narrow range ±0.1i"T in the vicinity of 
2. IK, where p n /p — ► 1. 

In conclusion, it should be emphasized that the ob- 
served temperature (actually, frequency) dependence of 
the snowball effective mass is interesting not only in it- 
self. It also turns out sensitive indicator revealing quali- 
tative difference in the velocity distribution field around 
the moving sphere in either viscous or ideal regime. This 
difference has been known for a long time. However, it 
did not attract much attention since slow decrease of the 
velocity field (~ 1/r) does not result in any divergence 
in the real part of the Stokes drag force. The situation 
is quite different in calculation of its imaginary part, and 
outlined above insight is one of important conclusions of 
this paper. 

Certainly, the hydrodynamic treatment of the M(T) 
dependence imposes some restrictions on explanation of 
the available experimental data obtained at low temper- 
atures (in the vicinity of several mK). However, the al- 
ternative kinetic language suitable for description of the 
ballistic regime confronts in substantial difficulties when 
the sphere effective mass is calculated even in the lami- 
nar limit. Details of this formalism will be reported else- 
where. 
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